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ABSTRACT 

One scenario proposed to explain the million degrees solar corona is a finely- 
stranded corona where each strand is heated by a rapid pulse. However, such 
fine structure has neither been resolved through direct imaging observations nor 
conclusively shown through indirect observations of extended superhot plasma. 
Recently it has been shown that the observed difference in appearance of cool and 
warm coronal loops (~ 1 MK, ~ 2 — 3 MK, respectively) - warm loops appearing 
"fuzzier" than cool loops - can be explained by models of loops composed of 
subarcsecond strands, which are impulsively heated up to ~ 10 MK. That work 
predicts that images of hot coronal loops (> 6 MK) should again show fine 
structure. Here we show that the predicted effect is indeed widely observed in an 
active region with the Solar Dynamics Observatory, thus supporting a scenario 
where impulsive heating of fine loop strands plays an important role in powering 
the active corona. 

Subject headings: Sun: corona — Sun: UV radiation 



1. Introduction 

The bright corona consists of magnetic loop-like tubes which confine the heated plasma. 
It has been proposed that the plasma is heated by rapid energy pulses, the so-called nanoflares, 
due to very localized reconnections of the magne tic field braided and twisted by the chaotic 
motio ns of the loop footpoints in the photosphere ( )Parkerlll988l ; ICarg i ll994i ICargill fc Klimchukl 



2004 ). Although impulsive events were predicted to be small and rapid, they are expected 
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to produce hot, variable emission at ~ 10 MK, even in the absence of flares. The ev i- 



dence for these si g natures has been nei t her direct nor conclusi ve so far ( Reale et al.ll2009bl 



McTiernanl 120091 ; ISchmelz et al.l l2009t ISylwester et al.l 120101 ) . The lack of evidence sug- 



gests that alternative mechanisms, i.e . the more gradual di ssipation of MHD Alfven waves 



(lHolrweglll984j ; iNakariakov et al. 



1999 



Of man fc Wane 



20081 ) , or even the direct involvement 



of the underlying chromosphere (IDe Pontieu et al.ll201ll ). are important (even critical) to the 
solution of the heating problem. 

Indeed, there are good rea sons why extensive nanoflaring activity has been elusive so 
far (lKlimchukl|2006l ; lRealdl2010l ). An important one is that, since momentum and energy are 
transported only along the magnetic field lines, a single pulse heats just one long and thin 
strand at a time. These strands are too thin to be spatially resolved by curr ent telescopes. 



and t he detection of hot plasma is difficult because of its small filling factor (IMartens et al. 
1985h . 



Recently, it has been shown that coronal loops substructured in a multitude of thin 
strands pulse-heated up to 10 MK can explain the e vidence of increasing fuzziness with 
temperature of emitting plasma (IGuarrasi et al.ll2010l ). This model also predicts that the 
fuzziness will decrease in bands sensitive to temperature > 3 MK. Here we show that this 
prediction is indeed confirmed by active region observations with the Solar Dynamics Ob- 
servatory. 

We suppose that, in the multi-stranded loops that populate an active region, the du- 
ration of each energy pul se is much smalle r than the plasma characteristic radiative and 
conductive cooling times (ISerio et al.l 1 19911 ). Once the heat pulse ends, the plasma cools 
down exponentially. If each loop strand is heated infrequently, i.e. the time between heating 
events is comparable or larger than the radiative cooling time, the plasma in the strand will 
be at high temperature only for a short time, and cooling for most of the time. Therefore, 
at a given time, only a few strands will be very hot, and many of them will be significantly 
cooler than their peak temperature. If nanoflares are scaled-down versions of coronal flares, 
we expect peak temperatures of about 10 MK. On the other hand, while the plasma is so 
hot, its emission is faint because of its low density: dense (and bright) plasma comes up from 
the chromosphere on time scales longer than the duration of the heat pulses. Therefore, the 
strands become bright when the plasma is already cooling and remain bright until the plasma 
has drained significantly. More specifically, it has been shown that a strand heated to 10 MK 
for ~ 1 min is at its brightest when it has cooled to ~ 3 MK, and that it remains around 



this temperature for a relatively long time ( IGuarrasi et al.ll2010l ). Thus, in this scenario, we 
expect to see few hot 10 MK strands, and many 3 MK strands in active regions. 
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2. Data analysis and modeling 



The Solar Dynamics Observatory (SDO) was launched in February 2010. The Atmo- 
spheric Imaging Assembly (AIA), with its 7 EUV narrowband channels, images the solar 
corona with high spatial resolution (~ 0.6arcsec/pixel), and high cadence in several spectral 
bands at the same time. We analyze SDO/AIA observations of AR11117 on 28 October 2010, 
from 2:00 UT. We consider a 500 x 500 pixels region in three channels: 17lA, 335A, 94A. 
The narrow passbands of the AIA channels contain bright spectral lines emitted by plasma 
at different temperatures. In particular, the channels centered at 94A and 335A contain 
a strong line of FeXVIII and Fe XVI, respectively, wh ich are emitted more efficiently by 
plasma around 8 MK and 3 MK ( lO'Dwyer et al.l 120101 ) . Therefore, these two channels are 
appropriate to test the prediction of the model. 




Temperature (K) 



Fig. 1. — Instrument response per unit emission measure and as a function of temperature, for the 
94A, 335A and 171 A SDO/AIA channels. 

We use level- 1.0 data, after standard processing of level-0 data (bad-pixel removal, de- 
spiking, flat-fielding). Data are obtained from a standard observing series with cadence of 
12 s in all channels, and exposure times of 2 s in the 171A channel and 2.9 s in the 335A 
and 94A channels. The channels have very different sensitivity to the solar coronal emis- 
sion (Fig. [I]), resulting in different signal-to- noise ratios. Therefore, for a more meaningful 
comparison of the active region morphology in the different channels, and to avoid possible 
spurious effects due to the different noise level, we summed images in the lower intensity 
channels: first we coaligned the images in each channel by using a standard cross-correlation 
routine (tr_get_disp.pro in IDL SolarSoftware package); then we added 30 consecutive images 
in the 94A channel, and 3 in the 335A channel (noise of ~ 3 — 4% in both channels in the 
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central region of Fig. [2]). Final images in all three channels were then co-aligned again. We 
note that during the 94A time interval (~ 7 min) the region variability is of a few percent 
at most, and that our conclusions on the observed morphological differences in the hotter 
vs. cooler emission are, if anything, underestimated by the temporal averaging we effectively 
apply by summing up images. 

Figure [2] (top three images) shows the internal part of the active region in the 94A, 
335A and 171A channels. The 171A channel contains a very strong FelX line with peak 
emissivity at ~ 0.8 MK. For the moment, we focus our attention to the top two images, and, 
in particular, to the core of the active region. 

In the 335A channel, the active region core is covered quite uniformly by a large number 
of bright arches. In the southern part (20 < Y < 55 arcsec), the arches coalesce to form 
a uniform bright band. In the northern part, we can identify three brighter loop groups in 
a background of more diffuse emission. A few very bright spots are visible at the center of 
the region. Overall the region has quite a diffuse appearance and individual loops cannot be 
clearly resolved. 

In the 94A channel, overall we see a very similar morphology and many bright structures 
are clearly cospatial with those observed in the 335A channel. The most striking difference 
from the image in the other channel is in the core itself: in the southern part, while in 335 A 
the arches are densely packed and uniform, in 94A they have greater contrast, i.e. we see 
an alternation of bright and fainter structures. In the northern part, we are even able to 
resolve very thin bright east-to-west bridges, in the same location where thicker arches are 
present in the softer channel. Overall, in the 94A channel, the loop systems appear sharper, 
the observed emission largely less "fuzzy" , and we can resolve thinner bright structures than 
i n the 335A channel. This is exactly the effect that we expected, and that was predicted 



flOuarrasi et al.ll2010h 



We have however to be cautious in one important point. Although the channel passbands 
are narrow, they include several spectral lines. In particular, the 94A channel includes 
another strong line (FeX) which peaks at ~ 1 MK. In general, we cannot be sure that the 
emission imaged by this channel comes only from hot plasma (> 6 MK). The 171 A image 
helps us in this respect, because it allows us to localize the bright cooler plasma (T ~ 1 MK), 
and to assess whether the 94A emission is due to hot or cool plasma. In the 171 A channel 
the active region shows quite a different morphology. Many structures are complementary 



to those observed in the other channels (jReale et al.l 120071 ). The core appears depleted 



of arch-like structures. Only few of them are visible, and they look quite different from 
those in both the other channels. The arch-like structures are instead replaced by bright 
"moss" . This moss is a well-known feature of this soft channel, already studied in Normal 



- 5 - 






20 40 60 80 100 

X (arcsec) 

Fig. 2. — Inner part of the active region AR11117 observed in 3 different channels (335A, 94A, 
171 A) of the SDO/AIA on 27 October 2010 around 02 UT. The color scales as the square root of the 
pixel counts. The ranges are 59-1567, 55-3157, 155-4681 DN for the three channels, respectively. 
The channels are most sensitive to plasma emitting at the labeled temperatures. The bottom panel 
shows the image obtained by subtracting the cool component scaled from the 17lA channel to 
the 94A image (same color scale as second panel, real range 0-2741 DN). Fig. [3] shows brightness 
profiles along the marked vertical and horizontal lines. 
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Incidence X-ray Telescope (NIXT) and Transition Region and Coronal Explorer (TRACE) 
observ ations and commonly explained as the bright warm footpoints of ( hot) high-pressure 
loops JPeres et al.lll994[ iFletcher fc de PontieJl999i iMartens et al.ll200oh . The 17lA image 



clearly indicates that much of the plasma confined in the filamented arches that we see in 
the 94A channel is not warm at 1 MK, and therefore it must be hot around 6-8 MK. As 
a further test, we estimated the emission measure of the cool plasma from the 17lA map, 
and used it to compute the expected contribution of the cool plasma in the 94A channel. 
We made the conservative assumptions that the cool plasma temperature is the one of the 
peak of the cool 94A response component (Fig. [TJ), i.e. logT = 6.05, a nd that the cool 94A 



response component is underestimated by a factor 5 (I Aschwandenl |2 1 ll ) . We subtracted this 
estimated cool 94A map from the 94A map shown in Fig. [2j The resulting image is also shown 
in Fig. [^bottom) . We note that the moss emission visible in the 94A channel is considerably 
reduced, and that the loop footpoints regions and the outbound loops are mostly removed. 
On the other hand, the emission of the loop bundle in the central region does not change 
significantly, except for a reduction in the moss region, and the fine structures in the upper 
half of the image are also almost unchanged. This confirms that the emission of the loop 
bundle and of the fine structures mostly comes from hot plasma. 

For a quantitative estimate of the different fuzziness of the hot (> 6 MK) and cooler 
(~ 3 MK) plasma, Figure [3^,b plots the pixel brightness, divided by a moving average, 
along the vertical and horizontal lines in Figure O The brightness of the 94A channel after 
subtracting the cool 1 MK component is also shown. We note that all the steep gradients are 
traced by more than one data point, i.e. the trends are coherent and are due to real structures 
not to noise. Along the horizontal line - that runs approximately along the magnetic tubes - 
the pixel brightness changes in a similar way and with a similar amplitude in both channels. 
Table [1] shows the root-mean-square (RMS) averages of the amplitude excursion, i.e. ~ 3% 
and ~ 5 — 6%, respectively. Along the vertical line - which runs across the field line direction 
- both brightnesses in the 94A channel are significantly more variable in space than in the 
335A channel, and the RMS average of the amplitude excursion becomes 10-13% vs 6%. 

To check on larger baselines, we derived the analogous values over all the rows (X) and 
columns (Y) of selected regions. For each region, Table [T] shows the RMS average values, 
and their standard deviations, in both channels. We selected three regions: one with as little 
moss as possible, the very center, and a larger core region. For the second region, we also 
report the values obtained in the 94A channel after subtracting the cool component. We 
find that the ^-values in the 94A channel are invariably much higher than all the others, 
thus confirming that the decrease of fuzziness in the 94A channel is highly significant. 

To add support to this evidence, we take the 2-D Fast Fourier Transforms (FFT) of the 
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Fig. 3. — Brightness profiles in the 94A before (red) and after (black dotted) subtracting the cool 
1 MK component, and in the 335A (green) channels along the horizontal (a) and vertical (b) lines 
in Fig. [2j The profiles are normalized to a moving average with a 10 pixels boxcar. A typical error 
bar is also shown, (c) Normalized 2-D Fourier Transforms of 335A and both 94A images in Fig. [2j 
94A transforms are both systematically higher than the 335A transform in the wavenumber range 
between the dashed lines. 
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images of Fig. [2] in the 94A channel before and after subtracting the cool component and 
in the 335A channel. Then we su m along circles of co nstant wavenumber, as it was done 



in the past for NIXT observations (IGomez et al.lll993al lbl). The power at zero wavenumber 



reflects the mean intensity of the image, so we normalize the images by the average intensity 
to make cross comparison of the power distributions more straightforward. In Figure [3fc the 
resulting Fourier transforms are power-laws. We see that both 94A FFTs are systematically 
higher than the 335A FFT in the wavenumber range ~ 10 — 70, that corresponds to a spatial 
range ~ 14 — 2 pixels. This range is in agreement with the cross-section size of the observed 
structures. All the lines converge at the highest wavenumbers, which are related to the 
presence of noise. This indicates a similar signal-to-noise ratio for both channels. In the 
end, the figure confirms more power at relatively high wavenumbers in the 94A channel both 
with and without subtracting the 1 MK component. 

For reference, we now sh ow some results of hy dro dynamic modeling of loops consisting 



of thin pulse-heated strands (IGuarrasi et al.l 120101 ) . Each pulse lasts 60 s and brings the 



plasma temporarily to ~ 10 MK. The plasma then cools off freely. Images comparable to 
the observed ones are obtained, for instance, by assuming that an instrument resolution 
element includes about 60 unresolved strands. Each strand is symmetric with respect to 
the loop apex. Figure H] shows the stranded-loop emission obtained by folding the density 
and temperature profiles from the hydrodynamic model with the 94A and 335A channels 
responses. Only the coronal part of the loop (upper 75% of the loop length) is shown. To 
emphasize the effect we are addressing, we normalize the emission in each pixel strip across 
the loop. In Fig. H] we can clearly see more white and black thin strips in the 94A image, 
while the other is more uniformly colored. 

The root mean square average of the amplitude excursion obtained from this model is 
17% for 94A channel, and 9% for 335 A channel, not too far from those measured in the 



Table 1. Fractional RMS amplitude excursions 



Region 


335A X 


a 


94A X 


(7 


335A Y 


a 


94A Y 


<7 


Single row/column 


0.03 




0.05 




0.06 




0.10 




Region little moss (size: 60 X 60 pixels, 


0.032 


0.017 


0.046 


0.021 


0.052 


0.011 


0.082 


0.017 


center: [X = 58", Y = 38"]) 


















Central core region (70 X 70, [68", 49"]) 


0.044 


0.018 


0.059 


0.021 


0.072 


0.016 


0.111 


0.036 


Central core region (94A-171A) 


0.044 


0.018 


0.068 


0.025 


0.072 


0.016 


0.130 


0.037 


Whole core region (100 X 100, [72", 48"]) 


0.053 


0.021 


0.065 


0.021 


0.063 


0.011 


0.099 


0.020 



Note. - 



The first row is computed along the lines marked in Fig. [2] a is the standard deviation of the rms excursions. 




Fig. 4. — Emission of a (straightened) pulse-heated multi-stranded loop synthesized from hydro- 
dynamic modeling in the (a) 94A and (b) 335A channels. The emission in each pixel strip is 
normalized across the loop. The color scale is linear between 0.6 and 1.4 (average 1): white is the 
maximum of emission, and black the minimum. Only the coronal part of the loop (upper 75% of 
the loop length) is shown. 



data. As explained by iGuarrasi et al.l (120101 ). the ultimate reason for this behavior is that 
overall each strand spends proportionally a long time, and with a high emission measure, at 
temperatures around 3 MK, i.e. just the temperature of maximum sensitivity of the 335A 
channel. Extrapolating this result to a whole active region, threaded with thousands of 
magnetic field lines, the loops will therefore be more populated by bright strands and appear 
more uniformly bright in this channel. On the contrary, they will appear less "filled" in a 
hotter channel, consistently with what we observe in the 94A channel. 

Similar results are obtained with pulses that heat the plasma to slightly higher or lower 
temperatures (~ 20 — 30%). Simulations and model no longer match the observations for 
weaker pulses leading to plasma maximum temperature ~ 4 MK and lower, where the cool 
~ 1 MK component would dominate the observed emission in the 94A channel. We can 
also exclude models with heat pulses lasting significantly less than 1 min ute: the main hot 



emitt ing ions would not have the time to reach ionization equilibrium ( IReale &: Orlando 



2008 ). therefore effectively suppressing the high temperature emission, visible in the 94A 
channel. As a warning, we remark that this model is tailored to describe a system of some 
loops, rather than an entire active region, and that the center of an active region may require 
a different (larger) number of strands per pixel. 

As a further support, we only mention here that the model predicts the emission in 
the 94A channel to vary in time with larger amplitude than in the 335A channel, and that 
this is confirmed by the observation: the pixel light curves obtained from the model are in 
good agreement with the observed ones. We defer a detailed report on time variability to 
forthcoming work. 
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3. Conclusions 



The new SDO active region data show evidence of extensive hot plasma, which is very 
finely structured, as expected in a scenario where storms of intense and rapid energy pulses 
are heating multi-stranded coronal loops. The strands are temporarily heated to flare tem- 
peratures. Our analysis is straightforward, and little depends on processing details, such 
as background subtraction. We remark that in the 94A channel we are not resolving the 
individual strands. We instead still see bundles of strands and the brightest ones are those 
where the fraction of very hot strands is relatively large. The reason why these very hot com- 
ponents are so difficult to detect is their small emission measure, and their small duty cycle 
with respect to most of their evolution time, spent mo stly in the subsequ e nt long cooling 
phase. This result supports previous debated analyses ( Reale et al. 2009bl lal; ISchmelz et al. 



2009 



McTiernanl 120091 ; ISylwester et al.ll2010f ). is consistent with many other pie ces of evi- 
dence pointing to dynamically heated loops, e.g. overdensity of ~ 1 MK loops (IKlimchuk 



20061 ; lRealell2010l ). and supports localized heat pulses, e.g., nanoflares, to play an important 



role in powering coronal active regions. 
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